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Abstract 

We show that efficient quantum computation is possible using a 
disordered Heisenberg spin-chain with 'always-on' couplings. Such 
disorder occurs naturally in nanofabricated systems. Considering a 
simple chain setup, we show that an arbitrary two-qubit gate can be 
implemented using just three relaxations of a controlled qubit, which 
amounts to switching the on-site energy terms at most twenty-one 
times. 
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Many possible physical implementations for performing quantum 
computation have been proposed (see Ch.8 in pQ for example). How- 
ever the vast majority of these schemes require that the coupling be- 
tween qubits (e.g. two-level atoms or quantum dots) can be turned on 
and off. In addition to the difficulty of doing this sufficiently quickly 
so as to avoid decoherence, and sufficiently accurately so as not to 
introduce unwanted errors, there is the fundamental problem that 
coupling terms in the Hamiltonian of multi-qubit systems are perma- 
nently 'on', i.e. they are time-independent. Building upon the work of 
Refs. jHEl, Benjamin and Bose [I] recently made a significant advance 
beyond this paradigm by introducing a simple method for performing 
quantum computation with a Heisenberg spin-chain where the spin- 
spin couplings are 'always on', i.e. constant in time. This scheme 
is clearly of great potential use - however, a quantitative study of 
its efficiency and possible generalization has not yet been reported. 
Futhermore the couplings between each spin-spin pair were assumed 
to be identical jlj 03 El - however, this situation cannot be engineered 
reliably in systems such as nanostructure arrays. 

In this paper, we generalize Benjamin and Bose's scheme to a 
Heisenberg spin chain with non-identical couplings between spins (i.e. 
qubits). Our scheme of time- independent, non-identical qubit-qubit 
couplings, will arise naturally in nanofabricated systems such as ar- 
rays of quantum dots, or chains of Ceo buckyball cages. It can also be 
engineered to arise for atoms in optical traps. Within a simple chain 
setup, we show that an arbitrary two-qubit gate can be implemented 
reliably with at most three relaxations of one controlled qubit. Such 
two-qubit gates are a crucial ingredient for universal quantum compu- 
tation. This number of relaxations amounts to switching the on-site 
energies at most twenty-one times. The major difficulty within this 
type of quantum computing scheme arises when the controlled qubits 
are relaxed. However, our work shows that the number of times the 
controlled qubits are relaxed can be minimized using a simple search- 
ing procedure. Our findings should provide considerable comfort for 
experimentalists preparing to build candidate quantum computer sys- 
tems using coupled nanostructures: Time, effort and resources do not 
have to be wasted on engineering identical qubit-qubit couplings. 

The Hamiltonian we consider is the following: 

h = E E j^f + E ■>.,..! ■ mi C 1 ) 

3 3 

where j labels the qubit (e.g. quantum dot, Cgo buckyball or atom) 
and the qubit-qubit couplings are represented by {Jjj+i}. For sim- 
plicity we consider the case where the J's in each pair (Jj-ij, Jjj+i) 
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do not differ by more than ten percent in magnitude - however there is 
no indication that similar results will not hold outside this range. Re- 
call that single-qubit operations can be implemented if the following 
sum of terms 



is added to the above Hamiltonian. It is assumed that the experi- 
mentalist has control over the temporal form of the on-site energies 
{Ej(t), Fj(t),Gj(t)}. For example, these terms can be controlled by 
applying a magnetic field (hence changing the Zeeman energy) or elec- 
tric field (hence introducing a Stark shift) to the qubit j [Jj in the 
appropriate direction. We will henceforth assume that single-qubit 
operations are possible (which is essentially the same assumption as 
adopted in Refs. [13 E]) and will concentrate on two-qubit operations. 
Hence we will discard the single qubit terms given in Equation [21 since 
they are not needed for the discussion of two-qubit operations. 

Consider three arbitrary qubits within this chain. We assume that 
qubits 1 and 3 are work qubits and that the qubits and 4 adjacent to 
them are fixed. Without loss of generality, we can incorporate these 
interactions into the on-site energy terms of qubits 1 and 3. Focusing 
on qubits 1 to 3, we note that the four subspaces spanned by {|000)}, 



{|001), |010), 1 100)}, {|110), |101), |011>} and {|111» are invariant un- 
der H. Within the basis {|001), |010), |100)}, the operation H is: 




(2) 





where 



A3 — —E\ + E-i + E3 — J\2 + J23 • 
Within the basis { 1 110) , 1 101) , |011)}, the operation H is: 



A% — Ei + Ei — £3 + J12 — J23 

A2 = Ei — E2 + E3 — J\2 — J23 



(4) 
(5) 
(6) 



/ B x 2J 23 \ 
V= 2J 23 B 2 2J 12 
\ 2J 12 £ 3 J 



(7) 



where 



£1 
B 3 



Ei — E2 — E3 — J12 + J 2 3 



—Ei — E2 + E% + J12 — 
—Ei + E2 — E3 — J12 — 



J 23 

J23 



(8) 

(9) 
(10) 
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We assume that qubit 2 is initialized to |0). To perform a two-qubit op- 
eration on qubits 1 and 3, qubit 2 has to be decoupled from them after 
the operation. Namely, [e~ lUT ]n, [e~ tUT hi, [e~ iUr ]23 and [e _4t/T ]3i are 
all zero, where r is the duration of the operation and [e~ ]ij denotes 
the (i, j)-entry of e~ tUr in matrix form. The same requirements apply 
to the matrix [e _l ^ r ]. Due to unitarity of the operation, this amounts 
to requiring that the magnitudes of [e~ lJjT }22 and [e~ lVr }22 equal one. 
If this is achieved, then the corresponding two-qubit operation will be: 



/ 



N = 





-iUrl 



ill 





-iUr] 
-iU; 



J 13 






) 



(11) 



where 



a 



-i(E 1 +E 2 -E 3 +J 12 +J 23 )r 



e 22 



(12) 
(13) 



Our aim now is to deduce a set of ({Ej}, r) values which will allow 
us to decouple the controlled qubit after time evolution r. Namely, we 
are searching for solutions that satisfy the two constraints: {e~ lUT }22 = 
1 and [e- iVr ) 2 2 = 1, in a four dimensional space. The use of various 
optimization techniques in gate-building in the classical domain has 
been explored in Ref. [Hj. Here, we employ the gradient descent 
method with random initial positions, which is common to all control 
optimization problems (see Ch. 9 in Ref. [S| for example). Our search 
algorithm goes as follows: 

1. Given a range L and density parameter to, we divide the set 
{Ej} into m values {— L + :0<A;<m — 1} and we divide 
r into m values {L + ^\ : < k < m - 1}. 

2. For each of the m 4 point combinations, we calculate the corre- 
sponding matrices U and V. These constitute the initial position 
of our gradient descent method to find the parameters ({Ej}, r), 
with objective function: 



Obj({^},T) = 



-iUri 



22 



+ 



-iVr-i 



22 



(14) 



We then proceed with the gradient descent method. Specifically, 
we first fix r and perform a descent with respect to {Ej}, then 
we fix the set {Ej} and perform a descent on the coordinate 
r. We repeat the process until a local minimum is attained - 
specifically, we repeat the process until the change in the ojec- 
tive function is less than 0.000005. If the value of the objective 
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function at that local minimum is higher than 1.99995, we form 
the matrix N as defined above and calculate its classifying angles 
as in Refs. [TO1EEI]. 

3. If the angle is of the form (0, 0, c) with uncertainty ±0.00257r, it 
is recorded. 

Step 1 partitions the search problem into many possible starting po- 
sitions, in order that different two-qubit operations can be generated. 
The search in Step 2 ensures that the combination ({Ej}, r) decouples 
the controlled qubit from the work qubits as much as possible. We note 
that we separate the descents with respect to the set {Ej} and to the 
coordinate r because their effects on the objective function are very 
different. The set {Ej} affects the eigenvalues and eigenvectors of the 
matrices U and V while r acts as a scalar multiplier of the eigenvalues. 
The results are shown in Figure^with L = 5 and m = 30. For the cou- 
pling constants ( J i2 , J23) G {(1, 0.9), (1, 0.95), (1, 1.05), (1, 1.1)}, we 
find that almost all values of the rotation angle in the interval [0,tt/2] 
become represented. The largest gap between two nearest-neighbor 
points in terms of the rotation angle is ~ 0.00457T, while the mean gap 
between two nearest neighboring points in terms of the rotation angle 
is less than 0.00057T. It takes about 8 hours of computing time for each 
J value on a standard laptop computer. Since the algorithm is highly 
parallel, the computing time will be cut by one half if two computers 
are on and so forth. The particular coupling constants shown in Figure 
^are chosen to illustrate the algorithm - further simulations indicate 
that similar results should hold for other sets of values. In short, our 
simulations demonstrate that by controlling the on-site energies in a 

Heisenberg spin-chain, one can generate two-qubit gates of the form 
■ 12 

e CJo " z °" z with arbitrary c with some tolerable errors. Specifically, our 
simulations suggest the following conjecture: 

Conjecture. For all non-zero Jx2 an d J23, for any c G 
[0,7r/2] and for all e > 0, there exists a set ({Ej},r) such 
that 

Ti 2 e- iHi{E ^ T -e ciaZ ^ Z <e (15) 

where 1^2 (.) denotes the operation of tracing out the second 
qubit, H({Ej}) is as shown in Equation and \M\ returns 
the absolute value of the entry with the maximum norm in 
the matrix M. 

We have not been able to formally prove or disprove this conjec- 
ture. However we recall that our algorithm corresponds to searching 
for four-dimensional vectors ({Ej},r) which satisfy two constraints: 
[e~ tUr ]22 = [e~ iVT ]22 = 1- Hence we speculate that this excess in 
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degrees of freedom will enable us to perform different two-qubit oper- 
ations. 

References |1U1 11H I12| showed that an arbitrary entangling gate Q 
can be written as Q = SiAS r , where 

^4 = gCiwrf of e vxia\a\ gC 3 iof of 

with Si,S r G SU(2) ® SU(2). Since 

R k e cia ^R[ = e aa ? a 2 (17) 

where R k = e^/ 4 )*** <g> e^ 4 )*** G SU(2) ® SU(2) and fc G {X, Y}, we 
see that an arbitrary two-qubit gate can be implemented by letting the 
controlled qubits relax three times at most. This amounts to switching 
the on-site energy terms twenty-one times at most. For example, since 
the CNOT can be represented as (0, 0, 7r/4) it can be implemented 
in the present scheme with just one relaxation of the controlled qubit. 
We note that our measure of efficiency relates directly to the difficulties 
arising in physical implementations, i.e. the difficulties in relaxing the 
controlled qubits. Hence, the fewer number of times the controlled 
qubit need to be relaxed, the more efficient the computing scheme 
is. On the other hand, since non-trivial entangling only occurs when 
the controlled qubit is relaxed, one can easily relate our measure of 
efficiency to the number of entangling operations needed, which is the 
usual measure of efficiency in quantum computing (see Refs. j!21 113j 
for example). In this regard, our method compares well with the most 
recent result obtained in Ref. ^3] since our results also indicate an 
upper bound of three relaxations/entangling operations to implement 
an arbitrary two-qubit gate. 

Figure [21 shows a systematic plot of the switching profiles which 
illustrate this result. We note that the error in decoupling qubit 2 from 
qubit 1 and 3 could be further reduced by running longer simulations, 
as shown in Figure El The error in coverage of the angle range can 
also be improved with more simulation. Its presence can also be seen 
as making the quantum computation more probabilistic and hence 
should not be viewed as being detrimental. 

In conclusion, we have provided a generalization of the quantum 
computing scheme introduced in Ref. [I] and have also discussed its 
efficiency. We have found the important practical and theoretical re- 
sult that identical qubit-qubit couplings are not necessary for universal 
quantum computation, thereby eliminating a fierce requirement on the 
precision of nanofabrication. We believe that the results of this work 
are relevant to a range of quantum computing implementations, and 
should stimulate experimentalists to explore less mature fabrication 



6 



technologies. In particular, the flexibility of the objective function 
in the gradient descent step should allow experimentalists to weigh 
parameters differently in order to reflect particular experimental con- 
ditions, or characteristics of fabrication. For example if the length of 
time for temporal evolution is a major constraint, perhaps because of 
decoherence processes, a penalty term of the form — r could be added 
to the objective function in order to skew the optimization. 

In more general terms, the method introduced here is not restricted 
to the Heisenberg interaction Hamiltonian. Specifically, given a set of 
fixed parameters and a set of controllable parameters which would 
be determined by physical requirements, one could search over the 
controllable parameters to try to come up with specific gates, as in 
the present paper. Therefore we envisage usefulness for this scheme 
when applied to more general Hamiltonian systems. 

C.F.L. thanks University College (Oxford) and NSERC (Canada) 
for financial support. N.F.J, thanks the LINK-DTI (UK) project. The 
authors are also very grateful to Simon Benjamin for useful comments. 
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Figure 1: (color online) Results obtained using the algorithm in the main 
text with J 12 = 1, L = 5 and m = 30. A corresponds to the mean separation 
between two neighboring data points in terms of their classifying angles (left 
axis), O corresponds to largest separation between two neighboring data 
points in terms of their classifying angles (left axis); x corresponds to the 
number of data points collected (right axis). 
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Figure 2: (color online) Profile showing the switching of {E{\ as a function of 
time, in order to realize a generic two-qubit gate. Red dotted line is E\\ green 
broken line is Ei\ blue solid line is E3 The sequence shown requires a total of 
twenty-one switchings of the on-site energies. An energy value represented 
by the word 'Passive' is assumed to be sufficiently large such that qubit 2 is 
fixed at state |0). 
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Figure 3: Plot of a typical evolution of the objective function versus the 
number of descents performed. Final value of the objective function in this 
simulation is greater than 1.999999 and the classifying angles produced are 
(0,0,0.415971-). 
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